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ABSTRACT 

We calculate radio light curves produced by the bow shock that is likely to form in front of the 
G2 cloud when it penetrates the accretion disk of Sgr A*. The shock acceleration of the radio- 
emitting electrons is captured self-consistently by means of first-principles particle-in-cell 
simulations. We show that the radio luminosity is expected to reach maximum in early 2013, 
roughly a month after the bow shock crosses the orbit pericenter. We estimate the peak radio 
flux at 1.4 GHz to be 1.4 — 22 Jy depending on the assumed orbit orientation and parameters. 
We show that the most promising frequencies for radio observations are in the 0.1 < v < 
1 GHz range, for which the bow shock emission will be much stronger than the intrinsic radio 
flux for all the models considered. 

Key words: accretion, accretion disks, black hole physics, relativity, acceleration of particles, 
radiation mechanisms: non-thermal 
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1 INTRODUCTION 

The center of our Galaxy is known to host a black hole with a 
moderate mass of M = 4.3 x 10 6 M o (e.g.. iGenzel et alj|2010h . 
This low-luminosity black hole, which is identified with the ra- 
dio source Sgr A*, accretes at a very low accretion rate (see, e.g., 
Quataert et al. 1999), has not been observed to produce relativistic 
outflows, and only sporadically shows flares in X-rays (BaganofF 
et al. 2001). The quiet and unimpressive life of Sgr A* is slated to 
change in the late summer of 2013 when a dense cloud of gas w ill 
approach it along a highly eccentric orbit dGillessen et al.l2012al lbT). 
with periapsis well inside the Bondi radius of the central black hole, 
Rb ~ 0.04pc » 2 x 10 5 Rq, where the gravitational radius is defined 
as R q = GM/c 2 = 2 x 10~ 7 p c. 

iGillessen et alj J2012ah discovered the gas cloud, called G2, 
during their extended observational campaign of stars in the vicin- 
ity of th e Galactic Center. The y estimated the cloud's mass to be 
< 3 M m dGillessen et alj|2012al) and determined the orbital param- 
eters of its center of mass. A semi-major axis of a = 666 ± 39mas, 
eccentricity of e = 0.9664 ± 0.0026, and epoch of periastron of 
'o.cm = 2013.69 ± 0.04, yield a periapsis distance tha t is only 
4400 ± 600R G from the black hole dGillessen etZll2012bl) . The in- 
teraction of the cloud with the accretion flow presents a unique op- 
portunity to probe the local conditions of gas at radii R ~ 10 3 5 ^g- 
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This corresponds to a few percent of the Bondi radius of Sgr A* 
dYuan et alT 2003) and is a largely unexplored region of the accre- 
tion flow around the black hole. 

Using theoretical models of the accretion flow around Sgr A*, 
iNaravan et al.l d2012al) recently predicted that a bow shock will 
form ahead of G2 for a few months around the time of closest ap- 
proach. This shock will likely accelerate particles and produce de- 
tectable levels of synchrotron radiation. The expected parameters 
of the shock correspond to a highly unusual regime for particle ac- 
celeration, one that has hardly been studied until now. Although the 
shock velocity is non-relativistic, at a temperature of T ~ 10 8 5 K 
the upstream electrons are nearly relativistic. Therefore, the elec- 
trons are more easily accelerated to high Lorentz factors and the 
injection problem often encountered in shock acceleration is miti- 
gated. In addition, the unshocked gas will be already magnetized, 
leading to conditions that could be particularly ideal for seeing syn- 
chrotron radiation from particles accelerated by G2's bow shock. 

In this paper we extend the work by Naravan et al. (2012a) 
in several ways. We explore in more detail the propagation of G2 
through the accretion flow around Sgr A* using a general relativis- 
tic magnetohydrodynamic (GRMHD) numerical simulation of the 
pre-impact accretion flow (Narayan et al. 2012b). We consider var- 
ious orientations of the orbit of G2 with respect to the accretion 
flow and utilize the simulation to determine the physical parame- 
ters of the pre-shock gas. We perform first-principles particle-in- 
cell calculations of particle acceleration in the shock for the range 
of physical parameters sampled along the orbit. We consider the 
net effect of particle acceleration along the cloud's trajectory using 
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two limiting approximations: one in which the moving cloud col- 
lects and holds on to all the relativistic electrons accelerated by the 
bow shock (the plowing case) and one in which electrons acceler- 
ated in the shock are left behind and radiate in the unshocked mag- 
netic field (the local case). We calculate the spectrum of accelerated 
electrons along G2's trajectory and the corresponding synchrotron 
radio light curves for these two possibilities. 

The structure of the paper is as follows. In Section [2] we dis- 
cuss the adopted model for the accretion flow around Sgr A*. In 
Section[3]we discuss the parametrization of the cloud orbit with re- 
spect to the disk equatorial plane. In Section|4]we show results of 
our study. We begin in Section |4T| by calculating profiles of phys- 
ical parameters of the accretion flow along the cloud orbit. In Sec- 
tion !4.2l we discuss the physics of electron acceleration and present 
the particle-in-cell simulations that we performed for relevant gas 
and shock parameters in order to provide estimates of the accelera- 
tion efficiency and the shape of the spectrum of non-thermal elec- 
trons. In Section |4~3l we present the predicted radio light curves and 
the spectrum of the radio emission from the bow shock and discuss 
their dependence on the orbit orientation. In Section [5] we com- 
ment on the uncertainties and optimal observing strategies of the 
G2 impact. Finally, we summarize our work in Section [6] 



2 MODEL OF SGR A* ACCRETION DISK 

Accretion on to Sgr A* is believed to occur via an advection- 
dominated ac c retion flow (ADA F: Narayan & Yi 1994 , 19951; 
Naray an et al. GUI | Yuan et al. 120031: I Narayan & McClintockl 
20081) . This mode of accretion is present whenever the mass ac- 
cretion rate is low, roughly < 1% of Eddington; hence, it oper- 
ates in the vast majority of galactic nuclei in the Universe. In the 
case of Sgr A*, observations in radio, infrared and X-rays, coupled 
with theoretical models, have pr ovided partial inf ormation on the 
parameters of the accretion flow. lYuan et all d2003l) used the num- 
ber density of electrons at the Bondi radius n e (R B ) « 130 cnT 3 
(Baganoff et al. 2003), along with a model of the radio and submil- 
limeter emission from thermal electrons in the inner accretion disk, 
to estimate the mass accretion rate on to the black hole. They found 
M BH s4x 10~ 8 M yr~' . Other spectral models, based on radiative 
transfer calculations using post-processed numerical simulations of 
global accretion disks, led to somewhat low er estimates, M BH ~ 
2 x 10 - 9 M G yr 1 dMoscibrodzka et alj|2009l) . while measurements 
of Faraday rotation in the magnetized accreting gas gave values 
somewhere in betwe en the two estimates: M BH ~ 10~ 8 M G yr~' 
jMarrone et ai1l2007l) . 

Although constraining, the results above do not specify the 
entire structure of the accretion flow. In particular, while the mod- 
els are fairly reliable near the black hole and near the Bondi ra- 
dius, they are poorly constrained at the intermediate radii probed by 
G2's orbit. To fill in this gap, studies rely on numerical calculations. 
GRMHD simulations provide a proper treatment of the innermost 
regions as well as of the flow on larger scalesQ. Recently, a number 
of groups have studied radia tively inefficient accret i on flows using 
such n umerical models (e.j 
2003L lAnninos et al] 12005 
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For the work reported here, we make use of one of the ADAF mod- 
els around a non-spinning black hole (spin parameter a, = 0) de- 
scribed in INaravan et al.l (1201 2bl) . The simulated flow corresponds 
to a Magnetically Arrested Disk (MAD, Narayan et al. 2003), i.e., 
a system in which the poloidal magnetic field is strong enough to 
partially "arrest" the accreting gas. The MAD state is believed to 
be a reasonable approximation of most ADAFs in the Universe fed 
by gas with initially non-zero poloidal magnetic flux. The most im- 
portant feature of the particular MAD simulation under considera- 
tion is that i t was run for an unprec edentedly long time. The results 
reported in INaravan et al. (2012b) corresponded to a duration of 
l0 s GM/c 3 , but since then the simulation has been run further, up to 
a time of 2x 10 5 GM/c 3 . As a result, the accretion flow has achieved 
steady state out to a large radius ~ 200Rq- This is important for the 
present work, as we explain below. The fact that the simulation 
assumes zero black hole spin, whereas the black hole i n Sgr A* 
probably has a moderate non-zero spin dBroderick et al.lEoilh . is 
not important since spin affects only the innermost regions of the 
accretion flow, whereas G2 interacts with gas that is much farther 
out. 

At periapsis, G2 is estimate d to be located at a r adius of 
R = 4400W G from the black hole dGillessen etZll2012bl) . There- 
fore, in order to predict synchrotron emission from G2's bow shock, 
we need to know the properties of the ambient accreting gas at this 
radius. However, no global GRMHD simulation can hope to obtain 
direct estimates at such large radii because of the prohibitive com- 
putational times required. Because of this, we extrapolate the phys- 
ical quantities we obtained from smaller radii. Such an extrapola- 
tion is acceptable, provided that the simulation has reached steady 
state out to a large enough radius where the flow shows self-similar 
behavior (at small radii, the influence of the black hole causes large 
deviations from self-similarity.) This is what makes the MAD sim- 
ulation described above particularly good for the present applica- 
tion. Having reached steady state out to radii R ~ 200^0^ the A° w 
properties in this simulation have almost certainly entered the self- 
similar regime (see Narayan et al. 2012b). Hence, we may safely 
extrapolate out to the radii of interest to us. In the present work, we 
choose R = 150/?g as the anchor radius from which we carry out 
the extrapolation. 

We present in Figure[T|the radial profiles of density (top-left), 
temperature (top-right), azimuthal velocity (bottom-left) and mag- 
netic to gas pressure ratio (bottom-right). We plot the MAD simula- 
tion results for R < 150^g (left of the vertical solid line) and show 
power-law extrapolations at larger radii, R > 150/?g- We indicate 
the radius corresponding to G2's periapsis with the vertical dash- 
dotted line. The various colors and line thicknesses correspond to 
different values of the polar angle 9. 

We adjusted the density normalization to match the measured 
density at the Bondi radius, w e (^B) = 130cirT 3 (Baganoff et al. 
2003), using a power-law extrapolation R [ . The corresponding 
mass accretion rate at the black hole is M BH = 5 x 10~ 8 M o /yr, 
which is consistent with the estimates given above. The slope we 
chose for the density extrapolation differs from the canonical slope 
for an ADAF disk R~ 3/2 . However, as emphasized by Yuan et al. 
(2012a, b), an R l profile provides a better description of nearly all 
global ADAF simulations to dat^| This is true of our MAD sim- 



2012b; McKinney et al. 2012; Tchekhovskoy & McKinney 2012). 



1 If the innermost regions are not of interest and the important t imescales 
are sh ort, a non-relativistic approach, such as the treatments of |Pen et alj 
120031) and Chan et al. (2009) are sufficient. 



2 |Pen et all {2003) found a R~ - 72 scaling for the density in a set of large- 
scale non- r elativi stic MHD simulations. However, subsequent work by 
IPang et afl feplll) concluded that an scaling is a better description of 
non-convective magnetized ADAFs. 
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ulation as well. Especially in the vicinity of the equatorial plane, 
we see that the simulation results at smaller radii match seamlessly 
to an R l extrapolated power-law at larger radii. As expected, the 
largest densities are found at the equatorial plane (8 = n/2) and the 
smallest near the polar axis (9 = 0). 

The profiles of gas temperature show equally good conver- 
gence. These again show an R' dependence, indicating that the 
gas temperature scales proportional to the virial temperature. At a 
given radius, the hottest gas is above the disk surface, near the po- 
lar axis. The simulation does not distinguish between the ion and 
electron temperatures, whereas detailed ADAF models generally 
allow for a two temperature gas with T t j= T e (Narayan & Yi 1995). 
In most models, the plasma becom es two-temperature o nly at radii 
well below 10 3 i?o (e.g., Fig. 2 in lNaravanetai1ll995n . At larger 
radii, Coulomb collisions between protons and electrons are suf- 
ficiently frequent to enforce a single temperature. For the model 
of the accretion flow employed in the present work, the Coulomb 
equilibration time between protons and electrons at pericenter of 
G2 is ~ 10 6 s, which is shorter than the orbital time by a factor of 
~ 50 (and of course, shorter than the radial accretion time by an 
even larger factor). 

The azimuthal velocity scales with radius according to the ex- 
pected Keplerian dependence Rr [l1 . At a given radius, the largest 
azimuthal velocity is seen at the equatorial plane. The further from 
the disk plane, the slower is the rotation. We assume the radial and 
vertical (in 0) velocities to be zero. 

The magnetic to gas pressure ratio shows the worst conver- 
gence among the four quantities shown. Only inside the bulk of 
the disk does it settle down to a relatively radius-independent value 
~ .1, which is a typica l value for most MHD accretion disks (but 
sec iGaburov et af]|2012l for a discussion of the impact of the ini- 
tial magnetic field configuration on the saturated value of the pres- 
sure ratio). Off the mid-plane, the gas is strongly magnetized and 
the pressure ratio does not settle down to a constant value even at 
R = 150^g- Thus, extrapolation under the assumption of a constant 
ratio is not well justified. However, this will have relatively little 
impact on the light curves when compared to other uncertainties 
(see Section[5]l. 

For convenience, in Appendix [A] we give analytic formulae 
which reasonably approximate the disk structure for R > 150Rq- 



3 ORBIT ORIENTATION 

The orbi tal parameters of the c enter of mass of G2 are well con- 
strained (Gilless en et al. I l2012bh . In this work we assume that the 
shock follows precisely this orbit, neglecting possible hydrodynam- 
ical interactions between the cloud and disk gas which might affect 
the cloud front. Modeling these hydrodynamic effects will require 
MHD simulations that include both the moving cloud and the or- 
biting accretion flow, which we will present in future work. 

The cloud orbital plane is inclined with respect to the ob- 
server's line-of-sight and is roughly aligned with the well-known 
ring of stars around Sgr A* (Bartko et al. 2009; Levin & Be- 
loborodov 2003). The orientation of the accretion disk around 
Sgr A* is, however, poorly constrained. Therefore, in this paper, 
we allow the full range of possible relative orientations between 
the cloud orbital plane and the disk plane. Given the uncertainties 
in our model (see Section|5}, it is unlikely, but not impossible, that 
the observed lightcurve will put constraints on the disk orientation. 

Because of the assumed axisymmetry of the accretion flow, 
the following two angles completely specify the geometry: the in- 




Figure 2. Three-dimensional visualisation of the orbit orientation. The or- 
bital plane is inclined to the disk plane by an angle i. The argument of 
periapsis for the orbits shown is oj = 0, n/4 and n/2 for green, blue and 
purple orbits, respectively. The arrows denote the cloud direction of motion 
and the disk rotation. 



clination angle, (', and the argument of periapsis, a>. Figure|2]shows 
a three-dimensional visualisation of three sample orbits whose or- 
bital parameters correspond to those of G2. The inclination i corre- 
sponds to the angle between the orbital plane of G2 and the equa- 
torial plane of the accretion disk. Cloud orbits that counter-rotate 
with respect to the disk correspond to i < n/2, while co-rotating 
orbits correspond to n/2 < i < n. The argument of periapsis co is 
the angle between the line connecting the black hole and orbit peri- 
center and the line of nodes where the two planes cross. The orbits 
shown in the figure correspond to to = (green), n/4 (blue) and 
n/2 (purple line). 

The cloud crosses the disk equatorial plane twice during each 
full orbit. The locations and times of these crossings depend on to, 
but not i. Figure [5] shows as a function of co the epochs when the 
two crossings occur for the cloud center of mass (CM). For to = 
(green orbit in Fig. |2}, one of the crossings takes place exactly 
at periastron (t = 2013.69), while the other occurs half an orbital 
period, i.e., ~ 100 yr later. For a> = n/4, the CM crossings take 
place at t = 2013.57 and t = 2015.493 For to = n/2, the disk 
equatorial plane crossings are symmetric with respect to periastron 
and occur at t = 2013.33 and 2014.06. 

The particular choices of i and to that we considered in this 
work are listed in Table Q] where the crossing times are given rela- 
tive to the epoch of pericenter to. 



3 Note that, for to = -n/4, the crossing times will be reflected (in time) 
with respect to the epoch of periastron. Later, we describe results corre- 
sponding to to = n/4. Disk properties along the orbit for to = n/4 can be 
easily converted to to = -n/4 by reflection. However, the radio light curves 
presented in Section R31 are not related so simply. 
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Figure 1. Numerical disk model and its extrapolation. Panels show logarithms of density (top left), gas temperature (top right), azimuthal velocity (bottom 
left), and magnetic to gas pressure ratio (bottom right panel). Vertical lines show the radius of extrapolation (R = 150Kg, solid), the periapsis of the cloud orbit 
(R = 440(M?G: dot-dashed) and the Bondi radius (R = 2 X 10 5 Rq, dotted). is the polar angle and = njl corresponds to the equatorial plane. 



Table 1. Orbit orientation parameters 

Inclination (i) 

7r/3 counter-rotating 

7r/4 counter-rotating 

n/6 counter-rotating 
7i - rr/3 co-rotating 
n - 7r/4 co-rotating 
tt-tt/6 co-rotating 



Argument Disk equatorial plane crossings 

of periapsis (to) relative to the epoch of pericenter to 

single crossing at t = to 

tt/4 two asymmetric crossings at t = to — 0. 12yr and to + 1.8yr 

tt/2 symmetric crossings at t = to - 0.37yr and to + ().37yr 



4 RESULTS 

4.1 Accretion gas properties along the orbit 

The radio synchrotron emission that we expect as a function of time 
as G2 orbits around Sgr A* depends on the properties of the ambient 



accretion disk at various points along the orbit. Key quantities are 
the density, temperature and magnetic field strength in the ambient 
gas, and the Mach number of the cloud relative to the gas. Using the 
extrapolated profiles discussed in Section|2] it is straightforward to 
compute all these quantities along the orbit as a function of time 
for a given choice of the cloud orbital parameters, i and to. 

The top-most panels in Figure [4] show the variation of den- 
sity along the orbit. The columns correspond to inclination angles 
i = 7r/3 (left, high inclination; but also i = n - n/3, since the scalar 
fields are identical when reflected), tt/4 (middle), and n/6 (right, 
low inclination). In Figure[4] colors denote orbits with different val- 
ues of the argument of periapsis <±>. The horizontal axis shows time 
relative to the epoch of periastron to. The green lines correspond to 
to = 0, where the cloud crosses the disk equatorial plane only once, 
at pericenter. Consequently, all these density curves have a single 
maximum at t = to- Blue curves correspond to to = n/4, where 
the cloud crosses the disk plane twice. For the high inclination case 
there are two local maxima. The first one corresponds roughly to 
the first crossing of the equatorial plane (t = to - 0.12), while the 
second occurs much earlier than the second crossing (~ t + 0.3 vs 
t + 1.8) because of the fact that the decrease of density with in- 
creasing radius dominates over the effect of the latitude away from 
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Figure 3. Crossing time of the disk equatorial plane (in years) by the cloud 
center of mass as a function of the argument of periapsis ai. The vertical 
axis is non-uniform and follows the logarithm of r - fo,cM, where fo,cM 
is the epoch of pericenter for the center of mass. Points and labels denote 
crossing times for at = —n/2, —n/4, 0, tt/4, and n/2. 



the midplane. The second peak is not visible for low inclinations 
because in these cases the cloud never goes very far from the mid- 
plane in between the two crossings. For the same reason, all the 
density curves in the last column (low inclination case) lie almost 
on top of each other. Indeed, in the limit i — > (which corresponds 
to the orbital plane coinciding with the disk plane), a> will have no 
impact on any of the curves. 

The second row of panels shows the variation of temperature 
along the cloud orbit. As discussed in Section [2] at a given radius 
the temperature is the lowest at the mid-plane and increases towards 
the polar axis. As a result, for the highly inclined orbit shown in the 
left panel, the green curve in which the orbit crosses the disk plane 
exactly at pericenter has twin peaks, a few months before and after 
pericenter. This is when the cloud comes closest to the polar axis 
and encounters the hottest gas. The opposite behavior is seen for 
the purple curve with two symmetric disk plane crossings. Here, 
there is a single peak in the temperature exactly at pericenter, when 
the cloud moves away from the mid-plane. As in the case of the 
density, for orbits with low inclinations, all the curves are similar. 

The third row of panels shows the magnetic to gas pressure ra- 
tio. Close to the mid-plane this parameter is approximately equal to 
0.1, while closer to the poles it takes on larger values. Qualitatively, 
the behavior is similar to that of the temperature, except that, being 
a ratio of pressures, there is no variation with radius. The shapes 
of the various curves in the different panels are then easy to un- 
derstand. The fourth row of panels shows the magnetic pressure 
^mag = B 2 /8tt. This quantity does vary with radius and behaves 
quite a lot like the temperature. Thus, the curves show striking sim- 
ilarity to those shown in the second row of panels. 

Finally, the bottom row of panels shows the cloud Mach num- 
ber with respect to the disk gas. Here two sets of curves are plotted 



for each choice of ( and co, thick curves for counter-rotating orbits 
(inclination = i) and thin curves for co-rotating orbits (inclination 
= 7r - (')• As expected, the Mach number is lower for co-rotating 
orbits. Also, the difference between counter- and co-rotating orbits 
increases for lower inclinations (the two curves coincide for an or- 
bit perpendicular to the disk plane, i = n/2). The Mach number 
for counter-rotating orbits lies in the range 1 .8 - 3.0 while for co- 
rotating orbits it is in the range 1.3 — 2.0. Most importantly, we 
find that the Mach number is always greater than unity in all cases. 
This means that a bow shock is expected to form for every orien- 
tation of the orbit and thus we always expect some level of excess 
synchrotron emission from shock- accelerated electrons. 



4.2 Electron acceleration in the bow shock of G2 

A variety of astrophysical evidence suggests that particles are 
accelerated efficiently via the Fermi and shock-drift acceleration 
mechanisms in shocks (e.g., Blandford & Eichler 1987). Typically, 
these processes give rise to a non-thermal power-law tail in the en- 
ergy distribution of the particles. The parameters of the bow shock 
of G2 correspond to an interesting regime that has not been well 
studied. The shock is non-relativistic, but the upstream electrons 
are quasi-relativistic (assuming thermal equilibrium between elec- 
trons and protons), with temperature varying between T — 10 8 " A" 
and T - 10 87 A" along the cloud orbit (see the second row of panels 
in Fig. [4}- As a result of the relatively hot temperature in the accre- 
tion flow, the shock has a modest Mach number, ranging between 
M ~ 1.5 and M ~ 3 (see the bottom-most row of panels in Fig. |4}. 

We studied the acceleration of electrons in the bow shock of 
G2 by means of two-dimensional first-principles numerical simu- 
lations, with the particle-in-cell code TRISTAN-MP (Spitkovsky 
2005). The simulation setup parallels closely the one emplo yed by 
Riquelme & Spitkovsky (2011) and iNaravan et al.l J2012ah . with 
the magnetic field lying initially in the simulation plane, oriented 
at an oblique angle with respect to the flow velocity. We fixed 
the angle between the upstream field and the shock velocity to 
be i]/ = 63°, and the ratio of magnetic to gas pressure in the 
accretion flow to be P m - ig /P m = 0.1. Below, we argue that the 
same acceleration physics operates across a wide range of mag- 
netic obliquities, and for ratios of magnetic to gas pressure as large 
as Pmag/^gas = 0.5 (see the third row of panels in Fig. |4). We ex- 
plored the dependence of our results on the pre-shock temperature 
over the range 10 s 1 K up to 10 8 7 K, as appropriate for the region 
closest to the cloud pericenter where most of the emission will be 
produced, and on the shock Mach number, which we varied be- 
tween 2 and 3.5. 

For computational convenience, we chose a reduced mass ra- 
tio m p /m e = 100, but we verified that our results remain the same 
for larger mass ratios, when all the physical quantities are scaled 
appropriately. Specifically, we ran simulations spanning the range 
m p /m c = 25-400, fixing the pre-shock electron temperature (equal 
to the proton temperature), the shock Mach number and the ratio of 
magnetic to gas pressure, and we measured the time in units of the 
inverse proton cyclotron frequency toT 1 . We also checked the con- 
vergence of our results with respect to the spatial resolutior0 and 
the number of computational particles per cell (up to 128 per cell). 

We find that at the shock, a fraction of the incoming elec- 
trons are reflected backward by the shock-compressed magnetic 



We resolved the electron skin depth c/ai ps with 10 computational cells, 
and the electron Debye length /!rj c = ^JkT/m e c 2 c/oj pl . with a few cells. 
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Figure 4. Profiles of (from top to bottom) density, temperature, magnetic to gas pressure ratio, magnetic pressure and Mach number along the cloud orbit for 
inclinations: i = n/3 (left), jr/4 (middle) and n/6 (right column) as a function of time (fo denotes the epoch of pericenter of the shock front). Thick and thin 
lines in the bottom-most set of panels correspond to counter- and co-rotating orbits, respectively. Colors denote the value of the argument of periapsis at and 
correspond to the orbits shown in Fig. [2] 
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field (Matsukiyo et al. 201 1) or by scattering off of electron whistler 
waves excited in the shock transition layer (Riquelme & Spitkovsky 
2011). For quasi-relativistic electron temperatures, the reflected 
electrons are fast enough to remain ahead of the shock, resist- 
ing advection downstream by the oblique pre-shock field. As the 
shock-reflected electrons gyrate around the shock, they get ener- 
gized by shock-drift acceleration (e.g., Begelman & Kirk 1990) and 
form a local non-thermal population, just upstream of the shock. 
The constraint that the reflected electrons be fast enough to out- 
run the shock along the oblique field, thus participating in the 
shock-drift acceleration, can be recast as an upper limit on the 
field obliquity tf/. If the characteristic electron thermal velocity is 
Vth,e = y/kT/m e and the shock velocity is v^, electron acceleration 
is possible only if if/ < if/ ait , where the critical obliquity angle i/r clit 
satisfies «th, c cos^ cr j t = u s h- So, the criterion for efficient electron 
acceleration can be rewritten as 



^ ;£ lAcrit - arccos \M yjmjm p j , (1) 

showing that, even for our largest Mach number (M = 3.5), elec- 
tron acceleration will be a common by-product of the shock evo- 
lution, being allowed for all angles i]/ < 85°, if we take the re- 
alistic mass ratio m p /m e = 1836. For our reduced mass ratio 
m p /m c = 100, we chose the obliquity angle iff = 63° so that the 
constraint t// < (fr cr j t required for efficient electron acceleration is 
satisfied for all the shock Mach numbers t hat we investigated. 

As demonstrated by iNaravan et al. | d2012al) for a pre-shock 
flow with temperature T = 10 9 K, the electrons participating in the 
shock-drift acceleration process are accelerated ahead of the shock 
into a power-law tail with a slope p = 2.2 and containing roughly 
5% of the incoming electrons (i.e., the so-called acceleration ef- 
ficiency is j] as 0.05)0 The upper energy cutoff of the pre-shock 
electron spectrum steadily increases and the slope becomes flatter 
over time, suggesting that the distribution will asymptote at late 
times to a power-law tail with p < 2.2 extending to very large elec- 
tron en ergies (we refer to the left panel of Fig. 2 in INaravan et al.l 
(2012a) for details on the temporal evolution of the electron spec- 
trum). The electron energization at early times is governed by the 
shock-drift mechanism, but the diffusive Fermi process is expected 
to dominate at late times. We find that the counter-streaming be- 
tween the incoming flow and the shock-reflected electrons trig- 
gers the Weibel filamentation instability ahead of the shock (Weibel 
1959, Medvedev & Loeb 1999), with the wavevector perpendicular 
to the pre-shock field0 By scattering off of the turbulence gener- 
ated by the Weibel instability, the electrons pre-accelerated by the 
shock-drift mechanism will be injected into the Fermi process. 

If the electron acceleration in our simulations is still primarily 
controlled by the shock-drift mechanism, the electron energy gain 
Aym e r during a time At will be 

Ay m e c 2 ~ qE v e // At , (2) 

where E is the strength of the background motional electric field 
(originating from the drift of the upstream frozen-in magnetic field 
toward the shock) and v e // is the component of the electron velocity 
along the electric field. By taking v e // ~ v^e (or a few times larger 



than «th,e), we can rewrite the expression above as 

Aym.r lm„ 

— ~Msin<A J— (o d At, (3) 
kl \ m c 

which can be used to interpret the dependence of our results on the 
pre-shock conditions, as we now describe. 

In Figure [5] we show how the electron spectrum ahead of the 
shock depends on the upstream temperature (left panel, with fixed 
Mach number M = 3) and on the shock Mach number (right panel, 
with fixed temperature T = 10 S 5 K). In the left panel, we present 
the electron spectrum ahead of the shock for four values of the flow 
temperature: T = 10 SA K (orange solid curve), T = 10 83 ^ (yellow 
solid curve), T = 10 S3 K (green solid curve) and T = 10 87 X" (pur- 
ple solid curve). The location of the thermal peak in the electron 
spectrum scales linearly with the pre-shock temperature. When we 
rescale our spectra to account for this effect (i.e., we shift them 
along the x-axis by (T/10 K)' 1 ), the resulting distributions nearly 
overlap (compare the green solid line with the three dashed curves). 
The electron power-law tail starts at y — 1 - £kT/m c c 2 , where 
£ =! 7.5 is nearly insensitive to the pre-shock temperature (vertical 
green dotted line), and the tail contains about 5% of the incoming 
electrons. The fact that electron spectra for different upstream tem- 
peratures are identical, modulo an overall shift in the energy scale, 
can be easily understood from Eqs. ([T) and (|3}. The former demon- 
strates that the condition for efficient injection into the acceleration 
process does not depend explicitly on the flow temperature, at fixed 
Mach number. In addition, Eq. ([3) shows that the temporal evolu- 
tion of the high-energy spectral cutoff, normalized to the thermal 
peak, is insensitive to the upstream temperature. 

The dependence of the electron spectrum on the shock Mach 
number is illustrated in the right panel of Figure [5] showing that 
the non-thermal tail is most pronounced in high-A4 shocks. Yet, 
some evidence for electron acceleration to non-thermal energies is 
present even in the spectrum of JVt = 2 shocks (black solid line), 
whose high-energy end clearly departs from a Maxwellian distri- 
bution (black dotted line). We point out that the main difference 
among shocks with different Mach number is not in the normal- 
ization of the non-thermal tail or in the location of its low-energy 
end, but in the extent of the power-law component. Regardless of 
the Mach number, the tail contains a fraction ~ 5% of the incoming 
electrons, and its low-energy end is at y - 1 a £;kT/m c c 2 , where 
£ =! 7.5 is nearly insensitive to the shock Mach number (vertical 
green dotted line). The fact that the acceleration efficiency is inde- 
pendent of the Mach number follows from Eq. |[T), since i// <; i/r cr i t 
for all of our choices of Also, the increase in the high-energy 
spectral cutoff with Mach number is in agreement with Eq. ((3), stat- 
ing that the acceleration rate scales linearly with M. It follows that, 
at a given time (in units of w^ 1 ), the spectrum of shocks with higher 
Mach number will have stretched to more extreme energies. 

The non-thermal tail in the electron spectrum of M = 3.5 
shocks (red solid line) can be fitted as a power law of slope p a 2.4 
(red dashed line). We expect the power-law tail to become flat- 
ter with time, as its high-energy cutoff increases, a nd the spectral 
slope p should approach the value p s 2.2 found by Narav anet al.l 
fe012al) in a flow with T = W 9 K. In shocks with M < 3 (black, 
blue and green curves), given the limited extent of the power-law 



5 With time, a power-law tail of similar normalization and slope as the pre- 
shock spectrum will appeal' in the post-shock electron distribution. 

6 We point out that the Weibel mode can only be captured by means 
of multi-dimensional simulations and was, therefore, absent in the one- 
dimensional experiments of Matsukiyo et al. (2011). 



7 In shocks with higher Mach number (and fixed ijj), the condition if/ < 
<Acrit w iU break down, fewer electrons will be reflected back at the shock 
(Matsukiyo et al. 201 1), so a smaller fraction of the incoming electrons will 
be injected into the shock-drift acceleration process. This will approach the 
limit of cold upstream plasmas studied by Riquelme & Spitkovsky (201 1). 
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Figure 5. Electron energy spectrum just upstream of the shock, normalized to the pre-shock electron density, for parameters relevant to the orbit of G2. (Left) 
Electron energy spectrum at co c it =10.8 for different upstream temperatures, but fixed Mach number M = 3. The spectra plotted with dashed lines have been 
shifted along the A-axis by (T/10 8 ,5 K) , to normalize the energy scale. (Right) Electron energy spectrum at oj cs t = 10.8 for different Mach numbers, but fixed 
electron temperature T = lO 85 *'. In both panels, the vertical green dotted line marks the location of the low-energy end of the electron power-law tail, whose 
slope p ^ 2.4 is indicated as a red dashed line in the right panel. In the right panel, the black dotted curve shows a Maxwellian distribution with T = 10 s 5 K. 



tail, it is still premature to draw firm conclusions regarding the ap- 
propriate value of the power-law slope at late times. Below, we as- 
sume that, regardless of M, the power-law tail in steady state will 
have a slope p = 2.4. As discussed above, this is likely to be a 
conservative upper limit, since the power-law index at late times is 
expecte d to approach the value p = 2.2 obtained bv lNaravan et al.l 
J2012ah . 

We have also tested that neither the acceleration efficiency nor 
the spectral slope depend significantly on the ratio of magnetic to 
gas pressure (we have tried P mg /P gas = 0.1, 0.3 and 0.5) or the 
obliquity angle tf/, provided that tf/ < tf/ a n. In fact, from Eqs. (T) 
and 10, we expect that the ratio P mag /P glm should not appreciably 
change the spectral shape, whereas the field obliquity (in the range 
4> S lAciit) should only affect the evolution of the high-energy spec- 
tral cutoff, as illustrated in Eq. l[3}- 

4.3 Radio lightcurves 

As we showed in the previous sections, the bow shock that will 
form ahead of the cloud G2 as it moves supersonically through 
the ambient disk gas is likely to accelerate, via the mechanism de- 
scribed in the previous subsection, the already hot electrons to rela- 
tivistic energies. The original thermal distribution of electrons will 
thus be modified and a high energy power law tail will develop. 
When these high energy electrons gyrate magnetic field present in 
the disk, they emit synchrotron radiation. As we will show below, 
under favorable circumstances, the emission from the shocked elec- 
trons can overwhelm the steady-state radio emission from Sgr A*. 

In our calculations we assume the following parameters to cal- 
culate the total number of accelerated electrons. We take the cloud 
cross section (the transverse area of the shock perpendicular to the 
relative velocity between the cloud and the accretion flow) to be 
A = nR 2 , withtf = 10 15 cmrl and use an electron acceleration effi- 

8 R = 10 15 cm is the value assumed by Nara yan et alj fe012al) and corre- 
sponds to 2/3 of the measured half- width at half maximum of the cloud 



ciency equal to n = 0.05. We take the low-energy end of the power- 
law tail at y - 1 = £(kT /m^c 2 ), where £ = 7.5 and T is the electron 
temperature, and the slope of the high-energy tail as p = 2.4. 

Whether or not the shocked electrons will remain close to the 
shock is unclear. We, therefore, consider two distinct limiting cases. 
In one scenario, which we call the plowing model, we assume that 
all relativistic electrons remain in the shocked gas and radiate in 
the instantaneous shock-compressed magnetic field. In the second 
scenario, referred to as the local model, we assume that the shocked 
electrons are left behind at their original position in the accretion 
flow while the cloud and the shock move on, and that the electrons 
radiate in the original uncompressed field. We believe that the real 
situation will be somewhere between these two extremes. 



4.3.1 Results from the plowing model 

In the plowing model, all the accelerated electrons remain in the 
vicinity of the shock. Because of the long synchrotron cooling 
timescale (> 100 yrs, see Narayan et al. 2012a), there is no loss 
of accelerated particles and hence the number of radiating rela- 
tivistic electrons at each energy increases monotonically with time. 
However, this does not imply a monotonically increasing luminos- 
ity since the synchrotron emissivity varies with the magnetic field 
strength. The field in the shocked gas is proportional to the field in 
the ambient accretion flow, which varies as a function of position 
along the orbit. As discussed in the previous section (see Fig. @), 
after an initial increase, the magnetic field strength drops as the 
cloud moves away from the central black hole. Therefore, while 
we expect an initial rapid increase in the synchrotron luminosity 
(since both the number of relativistic electrons and the magnetic 
field strength increase as the cloud approaches periapsis), this will 



iGillessen et alj2012j) . The choice of the effective cloud cross section and 
the impact of changes in cross section along the the orbit are discussed in 
Section[5] 
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be followed by a slower decay (the number of shocked electrons 
will continue to increase but the magnetic field will drop with time). 

To calculate the lightcurve, we pre-compute the positions and 
velocities of the cloud as a function of time for given orbit angles 
i and to. For each differential time interval between t and / + At, 
we calculate the energy distribution of electrons accelerated in that 
interval, 

dN . . (p- l)(£kT/m e c 2 )"- [ 



dy, 



: AvAtn e n ■ for 7 - 1 > £{kT jrn^c ), 

(7 - \y> 

(4) 

where n e is the local disk number density and v is the velocity of the 
cloud relative to the gas. The differential distributions are summed 
to obtain the total distribution of electrons at any given time accord- 
ing to 



dN ^ dN 



(5) 



where the sum goes over all the time steps up to the given time t. 

The spectral synchrotron power is calculated using the 
standard formula for a power-law distribution of electrons 
dRvbicki & LightmarJI 19791) . 



PM= JP + 19\JP_± 
m e c 2 (p+\) \4 12) \4 12 



2nm e cv 
3qB 



where 



dN 

C= — fcy = 2). 
d7 

The strength of the magnetic field B is given by 
(7 + l)M 2 



, (6) 



(7) 



B = 



:B,\, 



(8) 



(7- l)M 2 + 2 

which accounts for the amplification of the ambient disk magnetic 
field B disk by a compression factor that depends on the shock Mach 
number M and the adiabatic index y = 5/3. Note that only the 
component of the field perpendicular to the shock normal is am- 
plified and that, for the sake of simplicity, we ignore the longitudi- 
nal component here. Taking into account the distance to Sgr A* as 
D = 8.3kpc, we obtain the observed synchrotron flux as 

Py(t) 



Fy(t) 



(9) 



4;r» 2 ' 

Figure [6] shows a set of lightcurves calculated with the above 
procedure for v = 1.4 GHz (where self-absorption is negligible, see 
Section |4l4"l >. The rows correspond to different inclination angles, 
starting from i = n/3 (top), through i = n/4 (middle), to i = n/6 
(bottom), while the columns correspond, from left to right, to dif- 
ferent arguments of periapsis, a> = 0, n/4, and n/2. The thick and 
thin lines refer to counter- and co-rotating orbits, respectively. All 
the curves show the expected behavior, i.e., a relatively sharp in- 
crease when the cloud approaches pericenter, followed by a slower 
decay when the cloud and the shocked electrons move to regions 
with a weaker magnetic field. For all orbit inclinations, the radio 
emission is predicted to start rising a few months before the shock 
reaches pericenter. Co-rotating orbits produce significantly weaker 
emission due to the lower relative velocity of the cloud with respect 

9 The underlying assumption is that the acceleration time is shorter than Af, 
and that electrons accelerated at earlier times do not modify the shock struc- 
ture at later times (the shock structure is only determined by the properties 
of the pre-shock flow at that time). 



to the gas, which results in a smaller number of accelerated elec- 
trons as well as a weaker shock-compressed field (since the shock 
Mach number is smaller). The maximal flux for counter-rotating 
orbits is in the range 13 - 22 Jy, with the exact value depending on 
the orbit orientation. For co-rotating orbits, the range is 3 - 7 Jy. In 
most cases, the maximal emission is occurs roughly a month after 
the shock reaches periapsis. For some orbit orientations two peaks 
are predicted in the light curve. 

4.3.2 Results from the local model 

In the local model, electrons are assumed to move out of the 
shocked region and to remain in their original location in the ac- 
cretion flow, radiating in the original uncompressed magnetic field. 
We neglect adiabatic expansion losses of the accelerated electrons, 
since the expansion time is typically larger than the evolution 
timescale of the cloud along its orbit. In this approximation, each 
accelerated electron radiates under the same conditions through- 
out the entire period of the cloud encounter; hence, the procedure 
described in the previous section has to be modified. Equation [4] 
which describes the energy distribution of electrons accelerated in 
a given time step, still holds, but instead of summing up the num- 
ber of electrons and calculating the emission, here we sum up the 
radiation from electrons at each location that are accelerated in all 
the time steps up to a given time t, i.e., 

N(t) 

P y (f) = 2>*taa(0. (10) 

71 = 1 

where P v , local (0 is calculated via Equation (|6j but with 

(r,y = 2). (11) 



C= — 

d 7 local 



This last quantity corresponds to the number of accelerated elec- 
trons in a given time step at time t and B = B disk is the original 
uncompressed disk magnetic field. 

By construction, in the local model, the emission increases 
monotonically with time. However, because the magnetic field is 
not amplified, the total emission is lower than in the plowing model. 
Figure|7]shows a set of light curves that exhibit exactly such behav- 
ior. All of the curves start to increase a few months before the shock 
reaches pericenter but none of them reaches a flux that significantly 
exceeds 4 Jy. The difference in the magnitude of the maximum flux 
between the plowing and local models, which is roughly a factor 
of 5, is caused by two factors. First, the magnetic field is weaker 
in the local model than in the plowing model, leading to a corre- 
sponding decrease in the synchrotron emissivity. Second, because 
the accelerated electrons are not accumulated in the local model, 
fewer electrons reach the most favorable radius for emission. De- 
pending on the orbit orientation, the maximal flux at 1 .4 GHz lies 
in the 1 .4 - 4. 1 Jy range. 

4.4 Spectra 

In the preceding calculations, we assumed that the energy distri- 
bution of the shocked electrons is described by a single power 
law with slope p = 2.4. This results in a power-law synchrotron 
spectrum (F v oc v") with spectral index a = —(p - l)/2 (Eq.[6]l, 
which implies that the flux is highest at low frequencies. However, 
at very low frequencies, the emission is limited by synchrotron self- 
absorption, which was estimated by Narayan et al. (2012a) to affect 
frequencies v < 1GHz. In this section, we calculate synchrotron 
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a> = 7r/4 



a) = n/2 




Figure 6. Radio lightcurves at 1.4 GHz for different orbit orientations assuming the plowing model for acceleration (Section 14.3.11 . Rows correspond to 
inclination angles i = rr/3, n/4 and rr/6 (top to bottom). Columns correspond to the argument of periapsis to = 0, tt/4 and n/2 (left to right). Thick and thin 
lines are for counter- and co-rotating orbits, respectively. The vertical line corresponds to the bow shock epoch of periastron tg. 



self-absorption for the parameters adopted in this work and plot the 
expected spectral energy distribution. 

Synchrotron emission at a given frequency is self-absorbed 
when its total intensity, calculated according to Equation (6J, ex- 
ceeds the intensity of the local blackbody emission at the character- 
istic temperature corresponding to the energy of electrons respon- 
sible for emission at that frequency. The emission is then limited to 
the blackbody value. 

Synchrotron emission peaks at frequency v for e lectrons with 
Lorentz factor y given bv dRvbicki & LightmarJll979h . 



Amn.cv 



r 



3qB ' 

which corresponds to a temperature T equal to 

yra e c 2 



3/c B 



(12) 



(13) 



The spectral power emitted by a black body with this temperature 



from a surface area A is 

/\,bb = AnB v {T), 
where B Y (T) is the intensity of blackbody radiation, 
2hv 3 1 



BAT) ■■ 



Jiv/kT . 



1 



(14) 



(15) 



The critical frequency below which self-absorption limits the syn- 
chrotron emission is given by the condition 



Py, 



(16) 



where P v is given by Equation l|6} or Equation dlOt . For simplicity, 
we will not treat the transition regime in detail. Instead, we will 
assume that the total spectrum can be represented as a broken power 
law, with a slope equal to — (p - l)/2 in the high-frequency range 
and equal to 5/2 in the low-frequency segment that is self-absorbed 
below the critical frequency. 

In the case of the plowing model, we assume that the radiating 
electrons are contained within a sphere of radius R = 10 15 cm and 
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Figure 7. Similar to Figure|6]but for the local model (Section [4.3.21 . 



A = 4nR 2 . For the local model, we correct the emission for self- 
absorption independently at each cloud location assuming that the 
differential emitting area is AA = 2nRAr, where Ar corresponds to 
the distance covered by the cloud between time t and t + At. Thus, 
the emitting area is effectively the surface of a tube of radius R 
traced by the cloud moving along its orbit. 

Figure [8] shows our calculated spectra at t = to + 0.05 for the 
most and least favorable orbit orientations considered in this work. 
The left and right panels correspond to the top-right and bottom- 
left panels in Figures|6]and|7] respectively, i.e., to orbit orientations 
producing the most and the least amount of emission. The solid 
lines in Figure [8] correspond to the plowing model and the dashed 
lines to the local model. As before, thick and thin lines represent 
the counter- and co-rotating orbits, respectively. The flux levels ob- 
tained with the various models extend over a wide range, reflecting 
our lack of knowledge of the orbit orientation as well as the bow 
shock structure. 

The magenta data points in Figure [8] show the intrinsic radio 
emission of Sgr A* compiled from the literature ( Davies et al .ll976t 
iFalcke & Markoftood : Izhao et al]|2003l ; lMarrone et al.ll2008h . If 
bow shock-related synchrotron radiation is to be unequivocally de- 



tected, the emission must exceed the quiescent emission of Sgr A*. 
It is clear that, at frequencies above 100 GHz, none of the models 
predicts a detectable level of emission. The maximum emission in 
the plowing model occurs at ~ 1GHz and the flux here is well above 
the intrinsic emission of Sgr A* . However, at the same frequency, 
the local model predicts the bow shock emission to be roughly com- 
parable to the intrinsic flux. All our models predict detectable emis- 
sion at levels above the intrinsic emission of Sgr A* at frequencies 
in the range 0.1 < v < 1 GHz. 



4.5 When is t Q ? 

All of the light curves and spectra we calculated so far are ex- 
pressed as a function of relative time with respect to the time t 
when the bow shock reaches pericenter. The exact time correspond- 
ing to to is unfortunately uncertain. Since the bow shock is expected 
to form in front of the cloud, t is certainly earlier than the time ? ,cm 
when the cloud CM reaches pericenter. However, the exact location 
of the bow shock is poorly constrained. As a rough estimate, we can 
assume that the shock forms at the location corresponding to the 
"surface" of the original spherical cloud at t = 2000.0, i.e., at a dis- 
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Figure 8. Spectra of bow shock synchrotron emission for most (i = jt/3, ai = n/2, left) and least (i = n/6, lo = n/2, right panel) favorable orbit orientation (out 
of the orientations considered) calculated at f = to + 0.05. Solid and dashed lines correspond to the plowing and local acceleration models, respectively. Thick 
lines show counter-rotating orbits while thin lines denote co-rotating orbits. Magenta points show the quiescent radio emission from Sgr A* . The vertical line 
shows the frequency v = 1 .4 GHz for which the lightcurves have been calculated. 



tance R = 1.5 x 10 15 cm dGillessen et al.|[2012al) ahead of the cloud 
CM along the orbit. This particular location precedes CM by ~ 5 
months. Therefore, we estimate the time at which the bow shock 
reaches pericenter to be t as 2013.2, i.e., the bow shock should 
reach pericenter in March 2013. Naturally, this estimate is rather 
approximate since the internal structure of the cloud is unknown. 



5 DISCUSSION 
5.1 Uncertainties 

Apart from the uncertainty in f discussed above, the model adopted 
here for calculating radio light curves relies on a number of other 
simplifications and assumptions. The most crucial of these is the 
adopted cross section of the shock region, which determines the 
number of shocked electrons and, therefore, the intensity of emis- 
sion. We assumed that the effective area is constant in time, A = 
nR 2 , with R = 10 15 cm correspondin g to 2/3 of the ha lf- width at 
half-maximum as measured by iGillessen e t al. (2012a). However, 
given that the cloud density is much larger than the density of the 
surrounding gas, it is possible that it will effectively plow gas with 
even larger area, thus increasing the magnitude of the flux enhance- 
ment. 

A related uncertainty arises from the assumed shape of the 
cloud. G2 is not expected to keep its spherical shape throughout 
the passage: simulations show that the tidal gravitational field will 
first elongate it along the orbit, then squeeze it in the o rbital plane 
at pericenter dGillessen et al . 2012a; S aitoh et al.l2012l) , and finally 
let it decompress. However, the spatial extent of the cloud is large 
and the whole cloud will not become compressed at the same time. 
In other words, when the center of mass reaches the pericenter, only 
the middle part of the cloud is squeezed while its front has already 
decompressed and the tail lags behind. 

To quantify the effects of tidal compression, we calculated 
the trajectories of 500 test particles following the approach of 
IGillessen et alj d2012al) . We adopted the set of parameters given 
therein: a spherical cloud at t = 2000.0 with Gaussian distribu- 
tion of density with full-width at half-maximum (FWHM) equal 
3 x 10 15 cm and a velocity dispersion cr v = 130km/s. At four 
epochs (/ = 2000.0, 2013.0, 2013.69 and 2014.5), we projected 



each particle onto the plane perpendicular to its nearest section of 
the cloud orbit and plotted their distribution in Figure|9] Due to the 
large density of the cloud when compared to the disk, the whole 
region covered by test particles may be considered likely to plow 
the gas. The black circle corresponding to R = 10 15 cm is always 
well within the projected cloud cross section. Our value is, there- 
fore, rather conservative and the actual radio fluxes may be signifi- 
cantly stronger than those given in Section |4~3l if the effective area 
is indeed larger. In such a case, the critical frequency for the self- 
absorption will decrease, leading to higher flux at low frequencies, 
which may make the bow-shock emission easier to observe in that 
frequency range. 

In Section |4~2l we discussed the mechanism of electron accel- 
eration. Given the long computational time required to obtain the 
power-law index of the accelerated electrons, we used only a coarse 
sampling of physical parameters throughout each potential orbit of 
the cloud. However, the impact of this coarse gridding is not over- 
whelming. For example, taking a power-law index p = 2.6 instead 
of p = 2.4 drives the calculated fluxes down only by a factor of 
~ 2 and does not alter our conclusions regarding detectability. On 
the other hand, the value p = 2.4 adopted here is likely to be a 
conservative upper limit. As described in section 4.2, the electron 
non-thermal tail should evolve toward flatter slopes at later times, 
approaching the value p = 2.2 adopted by Narayan et al (2012a). 



Finally, we note that there are uncertainties arising from the 
model for the accretion disk. While the gas temperature should al- 
ways be close to virial, the disk density is poorly constrained by di- 
rect observations between the known values at the Bondi radius and 
the innermost regions where there is an estimate of the accretion 
rate (see also Psaltis 2012). We adopted a density profile p oc R' 1 in 
between these two limits. In principle, the density may be a more 
complicated function of radius, which could lead to a different den- 
sity along the orbit than predicted here. Given that the radio fluxes 
computed here are directly proportional to the local density of the 
accretion flow, they can easily be scaled for other models of the 
accretion flow. 
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Figure 9. Cloud profile on a plane perpendicular to its orbit calculated with 5000 test particles for (left to right) t = 2000.0, 2013.0, 2013.69 (epoch of 
periastron) and 2014.5. The black circle shows the region within R = 10 15 cm corresponding to the assumed cloud crossection. 



5.2 Suggested Radio Monitoring Observations 

The calculations presented in this paper indicate that the observed 
radio emission from Sgr A* should undergo a significant bright- 
ening around or shortly after the pericentric passage of the bow 
shock, i.e., in spring 2013, assuming / corresponds to March 2013 
(see Section l4~5b . In fact, given the ~ 22 mas displacement of the 
expected bow shock from the center of the accretion flow, where 
the long-term radio emission originates, the bow shock emission 
should appear as a new radio source that is, in principle, resolvable 
from the existing emission (Narayan et al. 2012a). In addition, we 
showed here that, under both the plowing and the local accelera- 
tion model assumptions, the spectrum of the bow shock emission 
should be different from that of quiescent emission in a number of 
ways. First, the emission is expected to peak at lower frequencies, 
around ~ 0.3 - 1 GHz. Second, the bow shock should produce de- 
tectable emission at v < 1 GHz, where the quiescent radio emission 
is undetectable. 

The Karl V. Jansky Very Large Array (JVLA) is ideally posi- 
tioned to detect and monitor the radio flux from the bow shock of 
G2 at frequencies around 1 GHz and above. Our "local" case with 
the least favorable orbit orientation produces the lowest level of flux 
enhancement. It nevertheless predicts a flux of 1.4 Jy at 1.4 GHz, 
which is comparable to the current quiescent flux from Sgr A* at 
that frequency (see Figure 8). The plowing case, on the other hand, 
produces a radio signal that can be an order of magnitude or more 
higher than the current flux at 1.4 GHz but is shorter lived, with 
the flux enhancement decaying on a timescale of the order of a 
year. Monitoring the Galactic Center with the JVLA in 2013, there- 
fore, should be able to not only detect the new radio emission from 
the bow shock but also to clearly distinguish the plowing and local 
models we considered here. 

Given that the spectrum of the bow shock emission may also 
peak below 1 GHz for certain orbit orientations and acceleration 
models, as we discussed in Section 4, monitoring the Galactic Cen- 
ter at frequencies ranging from < 100 MHz to 1 GHz provides an- 
other optimal opportunity to probe the shock physics and the condi- 
tions of the accretion flow around pericenter. Furthermore, because 
Sgr A* does not normally produce emission at these low frequen- 
cies, the detection of a bright radio source will deliver a definitive 
signature of the emission associated with G2. To this end, the Giant 
Metrewave Radio Telescope (GMRT), which operates at 5 frequen- 
cies between 50 MHz and 610 MHz, as well as at 1.4 GHz, is ide- 
ally placed and provides a strong complement to the JVLA. Given 
its latitude of 19°N, it is at a highly favorable location to observe 



the Galactic Center. In addition, GMRT has a distinct advantage in 
collecting area owing to 30 dishes with 45 m diameter each. 

We encourage frequent monitoring at available bands between 
50 MHz and 100 GHz during 2013 in order to capture the increase 
in the flux from Sgr A* , to measure its rise time, and to detect the 
possible decay a year after pericentre. 



6 SUMMARY 

The G2 cloud is expected to penetrate the accretion flow near 
Sgr A* with a supersonic velocity and shock the gas in front of 
it. We showed that the shock parameters and the pre-shock electron 
temperatures are conducive to accelerating electrons to relativis- 
tic energies. These electrons can emit synchrotron radiation in the 
magnetic field of the accreting gas. 

In this paper, we presented detailed predictions for radio light 
curves during the passage of the cloud. We adopted a state-of-the- 
art global numerical solution of a radiatively inefficient accretion 
flow as a base for our disk model and obtained physical disk param- 
eters at large distances by extrapolating the numerical solution from 
within its converged, steady-state region. The parameters describ- 
ing the electron acceleration process were determined by a set of 
particle-in-cell numerical simulations set up for the expected con- 
ditions as a function of the flow temperature and the shock Mach 
number. We considered a set of possible orientations of the orbit 
with respect to the accretion disk and two models for the parti- 
cle acceleration: the plowing case, where all the electrons are kept 
in the shocked region and radiate in the shock-amplified magnetic 
field, and the local model, where the accelerated electrons are left 
behind the cloud front and radiate in their local, unshocked mag- 
netic fields. 

In both models, assuming the bow shock leads the cloud CM 
by half a year, we found that the radio luminosity is expected to 
reach peak values in spring 2013. For the plowing scenario this 
maximum will be followed by a decay, while for the local scenario 
the luminosity will stay constant due to the long cooling time of 
electrons. Calculated maximum luminosities at v = 1.4 GHz span 
the range 1 .4 - 22 Jy, depending on the adopted orbit orientation 
and acceleration mode, and are likely to exceed the intrinsic lu- 
minosity of Sgr A* . At lower frequencies, the predicted emission 
will overwhelm the intrinsic luminosity for all the models con- 
sidered by roughly an order of magnitude. We, therefore, recom- 
mend an observational campaign in 2013 at frequencies in the range 
100 MHz - 1GHz. 
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APPENDIX A: APPROXIMATED MODEL FOR SGR A* 
ACCRETION FLOW 

For convenience, we give here analytical fits to the vertical structure 
of the disk at radius R = 150Rg, based on the GRMHD simulation 
discussed in Section [2] These profiles may be used to anchor ex- 
trapolations of disk properties to larger radii R > 150/?g- Table I All 
gives expressions for vertical profiles of density, temperature, az- 
imuthal velocity and magnetic to gas pressure ratio together with 
the proposed slopes for radial extrapolation. Radial and vertical ve- 
locities may be neglected. Figure |ATJ compares the numerical so- 
lution (red solid) with the fits (black dashed lines). The fits were 
obtained by least squares fitting to functions corresponding to ex- 
pressions in Table lAll in the region 0.15 < 6 < it — 0.15. The region 
closest to the polar axis is least reliable and was not taken into ac- 
count. 
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Figure Al. Vertical profiles of (top to bottom) density, gas temperature, 
azimuthal velocity, and magnetic to gas pressure ratio at R = 150Sq for the 
numerical solution (red solid) and the analytical fit given in Table lAll (black 
dashed line). 6 is the polar angle and 8 = n/2 corresponds to the equatorial 
plane. 
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